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We study the critical current I c dependence on applied magnetic field H for multifacet 
YBa2Cu30r_5-Au-Nb ramp-type zigzag Josephson junctions. For many experiments one would 
like to apply a homogeneous field in the junction plane. However, even tiny misalignments can 
cause drastic deviations from homogeneity. We show this explicitly by measuring and analyzing 
Ic vs. H for an 8 facet junction, forming an array of 4 x (0- it-) - segments. The ramp angle 
is 6 r — 8°. H is applied under different angles 6 relative to the substrate plane and different 
angles <j> relative to the in-plane orientation of the zigzags. We find that a homogeneous flux 
distribution is only achieved for an angle 8^ ~ 1° — 2° and that even a small misalignment 
~ 0.1° relative to Oh can cause a substantial inhomogeneity of the flux density inside the 
junction, drastically altering its I c vs. H interference pattern. We also show, that there is a 
dead angle 9^ relative to 8h of similar magnitude, where the average flux density completely vanishes. 



PACS numbers: 74.50.+r, 85.25.Cp, 74.78.Fk 
A. Introduction 

Large (in a geometrical sense) Josephson junctions 
(JJs) are studied intensively since decades, e.g. in 
the context of Josephson fiuxon physics. Interesting 
and important physics is related to junctions where 
some regions obey the usual Josephson relation, while 
other regions incorporate an additional phase jump 
of 7r, which can be viewed as a negative critical 
supercurrent density j c < 0. These 0-7r junctions 
can be fabricated in different ways, e.g. by connect- 
ing a d-wave superconductor such as YBa 2 Cu307_,5 
(YBCO)i-i or Nd 2 -xGe x Gu0 4 -5 (NCCO^ to a 
conventional superconductor like Pb or Nb, or by 
using biepitaxial grain boundaries in YBCO^— , or 
superconductor-ferromagnet-superconductor (SFS) or 
superconductor-insulator- ferromagnet-superconductor 
(SIFS) junctionals. In the YBCO-Nb or NCCO-Nb 
structures the Josephson junction is often of the ramp 
type and the barrier forms a zigzag line parallel to 
the crystallographic a and b axes£~— , see Fig. [T] This 
type of junction has been important for determining 
the symmetry of the order parameter of the cuprate 
superconductors^^. 

In addition, several properties make zigzag junctions 
very interesting for ongoing studies. In particular, half- 
integer vortices (scmifluxons) can form spontaneously 
at the corners of the zigzag line and, thus, quasi-one- 
dimensional vortex crystals can be realized 3 -^. Also, un- 
der special conditions the zigzag junctions can be used 
to create <p o 16]17 or ip junctions, with an arbitrary value 
of the ground state phase ip^r— . Such junctions, if long 
compared to the Josephson penetration depth Xj, can 



carry mobile fractional vortices (splintered vortices) hav- 
ing many unusual propertie a 19 i 22 ~ — . Apart from these 
research areas, aiming mostly at long junctions, ramp 
junctions are also interesting for superconducting elec- 
tronics, e.g. in the context of self-biased RSFQ circuits 2 ^ 
or in the context of superconducting quantum interfer- 
ence filter o 27 ' 28 . In all cases a good understanding of the 
implications of the zigzag ramp geometry is required. 

Studies of Josephson junctions in many cases require 
the application of a magnetic field H which, in theoret- 
ical studies, is usually considered to be oriented "paral- 
lel" to the junction plane, leading to a homogeneous flux 
density ^.qH in the absence of self fields generated by 
the Josephson currents. Most Josephson junctions have 
a simple geometry where the barrier layer and also the 
superconducting layers are oriented parallel to the sub- 
strate plane. Then, "parallel" simply means parallel to 
the substrate plane (0 — 0°, using the coordinates defined 
in Fig.[T](a)). For ramp junctions the junction plane (bar- 
rier layer) is tilted relative to the substrate plane by the 
ramp angle r (0 r — 8° for the junction we study here), 
while the superconducting layers are partially parallel to 
the substrate plane and partially bent along the ramp, 
cf. Fig.Q](b). "Parallel" is thus not defined well. As it 
will be discussed in detail in Sec. C, there is an angle 
Oh, with < 0^ < r , where the applied field creates an 
almost homogeneous and uncompressed flux density in 
the junction. A magnetic field applied under this angle 
should be referred to as parallel. 

Experimental studies on the zigzag junctions have usu- 
ally been performed with the magnetic field applied per- 
pendicular to the substrate piano 2 - - — . One reason for this 
was the finding that, when aligning the field roughly par- 
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FIG. 1: (Color online). Sketch of ramp-zigzag junction: (a) 
whole junction with 8 facets, (b) cross section of one facet. 
Penetration and focusing effect of applied magnetic field along 
the ramp are indicated by black lines. Coordinates, and ramp 
angle 9 r , as well as angles 8 and <j> of applied field H are 
indicated. 

allel to the barrier layer, apart from a field scaling factor 
due to flux focusing, no essential difference to the per- 
pendicular field orientation was observed. Theoretical 
interference patterns I c vs. H, calculated under the as- 
sumption that the applied field causes a homogeneous 
flux density B in the junction plane, and the experimen- 
tal patterns agreed only qualitatively: the critical cur- 
rent I c was maximum when the flux per (0 - tt) - segment 
roughly equalled one flux quantum 2 - - — . However, in al- 
most any other respect experimental and theoretical I c 
vs. H curves were not even similar. 

For junction geometries where the barrier layer is ori- 
ented either parallel or perpendicular to the substrate 
plane it was shown previously that for a field orientation 
perpendicular to the substrate plane (9 = 90°) the flux 
density in the junction barrier becomes inhomogeneous 
in the case of a homogeneous applied fielder— . This 
leads, e.g., to a striking difference of I c vs. H for SIFS 
multifacet - 7r junctions (SIFS-MJ Js) , when measured in 
fields applied parallel and perpendicular to the substrate 
plan e 32 ' 34 . Moreover, at least for junctions with a barrier 
layer oriented parallel to the substrate plane (for which 
0h = 0) and H applied under an arbitrary angle 9 there 
is a "dead angle" 9d where the magnetic flux caused by 
the parallel (9 = 0°) and perpendicular (6 = 90°) com- 
ponents of H cancel, leading to a critical current which 
almost docs not modulate with i J 30 ' 32 i 35 . The dead an- 



gle 9d can be very close to zero, making proper junction 
alignment parallel to H very difficult, if not impossible. 

Obviously, ramp zigzag junctions arc considerably 
more complex than conventional overlap junctions or 
SIFS-MJJs and need separate consideration. A system- 
atic study under oblique fields seems necessary, having in 
mind that this type of junction is very useful for many fu- 
ture investigations. In the present study we have chosen 
a YBa 2 Cu 3 07-Au-Nb junction with 8 facets (4 x (0 - 7r -) - 
segments), each facet being 10 /im long. Wc have inves- 
tigated its I C (H) dependence as a function of 9 and also 
the in-plane angle (f>, where <j> — 0° corresponds to H ap- 
plied along the facets oriented in y direction (cf. Fig. [T]) . 
Below we show that the field component perpendicular 
to 9h leads to a periodically modulated flux density pro- 
file. The average value of the flux density caused by this 
component is enormously compressed, by a factor pa 100 
compared to B = [IqH . As a consequence, the dead angle 
9* d , measured relative to 9h, is very small, 9* d sa —0.38° 
for = 45° and 9* d sa -0.25° for = 0°. To achieve 
a more or less homogeneous flux density B, the field H 
must be aligned better than some 0.1° relative to 9h and, 
to achieve 1 <I>o or more per (0- n) - segment, [i H values 
of more than 3 mT are required. These conditions require 
quite dedicated experimental setups. Alternatively, real- 
istic theoretical analyses should be based on the case of 
9 sa 90°. 



B. Samples and Measurement Techniques 

The YBCO layer was grown by pulsed laser deposition 
(PLD) on a [001]-oriented SrTi0 3 (STO) single crystal 
substrate. The substrate temperature T s during the de- 
position of the 65 nm thick YBCO thin films was 770°C 
and the oxygen pressure Po 2 was 20 Pa. A 60 nm thick 
STO insulation layer on top of YBCO was also deposited 
by PLD at T s = 760°C and P 02 = 9 Pa. The targets were 
ablated using a KrF excimer laser at a repetition rate of 
2 Hz. After thin film deposition, the chamber was vented 
with oxygen up to 0.5 bar and the sample was cooled 
down to room temperature, with an annealing step at 
450°C for 30min. The zero-resistance transition tem- 
peratures T c of the YBCO films were between 88 and 
90 K. The YBCO/STO bilayer was patterned using op- 
tical lithography and Ar-ion milling under an angle of 
30° between surface normal and ion beam. To ensure 
a constant ramp angle for all junction orientations the 
sample was rotated about the axis normal to the surface 
during milling. With these parameters a ramp angle of 
9 r = 8° was obtained. After removal of the photoresist 
the sample surface was cleaned in a soft Ar rf-plasma, 
in situ. Before the next deposition steps, an interlayer 
of 4 unit cells of YBCO was grown by PLD with the 
same deposition and annealing parameters as before to 
provide recrystallization 3 ^. According to Ref. [36[, the 
thin YBCO interlayer is expected to become supercon- 
ducting only on the YBCO ramp area, but not on the 
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FIG. 2: (Color online). Current voltage characteristic of the 
8-facet ramp zigzag JJ at T — 4.2 K for 3 different values of 
magnetic field, applied at 9 = 129°: 14.3 /xT (main I c max- 
imum), 0//T and 37 /xT (7 C minimum). Inset shows current 
voltage characteristic at expanded scales. 
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FIG. 3: (Color online) Images of the 8-facet YBCO-Nb zigzag 
JJ: (a) SEM surface image. The d a ,2_j / 2-wave order param- 
eter of the YBCO layer is indicated by the cloverleaf struc- 
ture. White and black lobes are orientated along the crys- 
tallographic a and b axes and indicate the sign change of the 
order parameter. The dotted rectangle surrounds the ramp 
area of one facet, (b) LTSEM SV image taken at T = 8.0 K, I 
= 2.3 fj,A and fioH — — 0.34 /iT (central I c minimum at small 
offset field for = 90°). 



STO substrate and insulation layer. With electron beam 
evaporation a 9 nm Au barrier was deposited, followed 
by a sputtered Nb layer of thickness dm, = 100 nm. The 
plasma-cleaning and the last three depositions steps were 
done without breaking the vacuum. The Nb and Au layer 
were patterned by a final photo lithography and Ar-ion 
milling step. The resulting J J is sketched in Fig. [T] 

The samples were measured at T = 4.2 K in a magnet- 
ically and electrically shielded cryostat. The mounted 
sample typically had a misalignment of the applied mag- 
netic field relative to the substrate plane 9 s below 1°. 
An external field fioH of up to 3.5 mT could be applied 
and continuously rotated with respect to 9 via two per- 
pendicular coils operated in linear combination. The in- 
plane angle <f> was varied by mounting the sample with 
proper orientation relative to the in-plane field axis, re- 
sulting in a misalignment error (j) Q g ~ 2 — 3°. 

Below we discuss results from our most homogeneous 
sample. Fig. [2] shows current voltage characteristics, 
measured at three values of applied field. The current 
voltage characteristics were hysteretic for critical cur- 
rents > 2 /xA, with a junction resistance of R rj 6.6 0. 
To measure I c the bias current I was ramped up at fixed 
magnetic field until the junction switched to its resistive 
state. A voltage criterion V CI — 2 /xV was used to deter- 
mine I c , leading to a maximum over-estimate of I c by 
V CI /R ~ 0.3 /xA in the nonhysteretic regime. In the hys- 
teretic regime I c is underestimated by some 0.1 /x A due 
to premature thermal activation. 

Fig. [3] (a) shows a scanning electron microscopy (SEM) 
surface image of the YBCO-Nb zigzag junction. The dot- 
ted rectangle surrounds the ramp area of one facet, which 
is visible as a faint grey line. The top Nb electrode over- 
laps the YBCO ramp area by extra 3 /im (idle region). 



The ramp areas of all facets form a zigzag line parallel to 
the a, b axes of the YBCO film. 

To investigate the homogeneity of the critical cur- 
rent densities of individual facets j*, we have imaged 
the current distribution of the entire zigzag junction at 
zero magnetic flux using low temperature scanning elec- 
tron microscopy (LTSEM) . Details of the method can be 
found in Ref. In brief, the pulsed focused electron 
beam, which is scanned across the sample, causes local 
heating ST < 1 K on a lateral length scale of « 1-3 /jm, 
which determines the spatial resolution of this imaging 
technique. The measured integral quantity is the voltage 
V across the junction, which is biased slightly above I c . 
ST results in a local reduction of | j c (T) | . The correspond- 
ing change of the overall I c of the zigzag junction causes 
a slight change of SV(x,y), which depends on the beam 
position (x, y) on the sample surface. For a bias current 
slightly above I c and B = 0, SV(x,y) oc —j c (x,y). This 
imaging technique requires non-hysteretic current volt- 
age characteristics, which for our junction is not the case 
for T = 4.2 K, cf. Fig. H Therefore, the SV(x,y) image 
of the YBCO-Nb zigzag junction in Fig. [3] (b) has been 
taken at T = 8.0 K. The junction has been biased at I 
= 2.3 fx A. A field fioH = — 0.34/xT was applied to com- 
pensate the residual field in the LTSEM setup. The SV 
image clearly shows the alternating sign of supercurrent 
flow across neighboring facets. On the scale of the spa- 
tial resolution of our imaging technique neither defects 
nor j c asymmetries between facets are visible. Further- 
more, the critical current densities of the facets seems to 
be quite homogeneous. 
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FIG. 4: (Color online). Sketch of a single facet J J oriented 
along the x axis, as considered in Sec. C: (a) Cross section 
parallel to the (y, z) plane. Penetration and focusing effect 
along the ramp for H applied under 9 = 90° are indicated 
by black lines, (b) Cross section in the (x,z) plane at y — 
(W/2) cos 9 r with (c) expected flux density profile B(x) in the 
junction barrier. In the (y,z) plane B is oriented along 6 r . 



area and the resulting flux density is essentially oriented 
parallel to the ramp (9 = 9 r ). Further, the screening 
currents preventing magnetic field lines from crossing the 
Nb layer lead to a strongly enhanced flux density in the 
junction barrier, cf. Fig. [4] (a). Along x this effect is 
strongest at the center of the junction (x = L/2). It 
disappears at the edge of the Nb film, cf. Fig. @] (b, c). 

To account for these effects we may represent the re- 
sulting flux density profile inside the junction barrier, 
oriented along 9 — 9 r , as 

Bgo,f(x) = fi Hf 90 jN 90J [a 90 j + P 90 ,/(x)] (1) 

where the field compression factor /go,/ 3> 1. The nor- 
malization constant Ngoj is chosen such that the spa- 
tial average of Ngojfaaoj + P90,/Oe)) equals 1. The in- 
dex '90' indicates perpendicular direction of the applied 
field (9 = 90°) and the subscript 'f indicates the sin- 
gle facet scenario. The spatially dependent field profile 
Pgoj(x) > 0, to be specified later, is maximum in the 
center of the junction and is assumed to reach zero at its 
edges. Without idle regions Pgo,/ would drop to almost 
zero at the junction edges, i.e. ago,/ ~ 0. The effect of 
the idle regions is to make Bqqj nonzero here. This is 
parametrized by ago,/. 

Using the spatial average (Pgo,/) of Pgo,/(x), Pgo,/ can 
be regrouped as 



C. General Considerations 

1. Single-Facet Ramp Junction 

Before we address I C (H) of zigzag JJs let us first theo- 
retically consider a JJ with a single facet oriented parallel 
to the x axis, (see Fig.0]). The J J area extends from x = 
to a length L along the x direction. The bottom YBCO 
electrode grows in thickness along y, reaching its max- 
imum thickness dy^co = 65 nm at y — dy'Ico / tan(# r ) 
= 462 nm, which is the projection of the junction width 
W = dYBCo/ sin (^0 = 467 nm to the y axis. 

Let us first estimate how the magnetic flux density 
distributes inside the junction when H is applied per- 
pendicular to the substrate plane (9 = 90°). The situ- 
ation is sketched qualitatively in Fig. 2] We first note 
that dyBCO i s we ^ below the YBCO in-plane magnetic 
penetration depth Aybco ~ 150 nm. Near the edge of 
the YBCO film located at y = magnetic field lines can 
cross the YBCO film on the scale of the Pearl length^— 
A = 2A Y bco/^ybco ^ 700 nm. This implies that over 
the whole junction area field lines can penetrate the 
YBCO film freely for any angle 9. 

By contrast, the thickness of the Nb film is comparable 
to the respective penetration depth (ANb ~ 80 nm) and 
approximately constant over the ramp area. Some field 
lines will cross the Nb film at its left edge (cf. Fig. [4] 
(a)) on a length scale of ANb- However, due to the idle 
region this edge is far away from the junction area. Thus, 
no field line will cross the Nb film inside the junction 



Bgojix) = /jioHfgoj 1 + Pgo,/ (a:) 



(2) 



where Pgo./fc) = N 90J (P 90J (x) - (P90,/)) has vanishing 
average. 

The case of "parallel" fields requires some discussion. 
For theoretical considerations "parallel" should refer to 
an angle (relative to the substrate plane) where the flux 
density in the junction is homogeneous. We denote this 
direction as 9h- In the absence of idle regions a natural 
choice would be the ramp plane, i.e., 9} x — 9 r . Screening 
currents in the idle regions, however, will deform the flux 
density similar to the case of perpendicular fields, leading 
to an ambiguity of what "parallel" actually means. 

Let us consider a field H applied at 9 = 0°, i.e. parallel 
to the substrate plane. Screening currents in the vicinity 
of the ramp cause B to be (almost) parallel to the Nb 
film, i.e. tilted by 9 r within the area of the YBCO ramp. 
These screening currents also cause a slight field com- 
pression towards the center of the ramp. For the case 
of 9 = 0° we thus expect a flux density profile for the 
component along 9 r of the form 



Bo, f {x) = tioHfojNoj [a , f + P ,/(x)] 



(3) 



where the subscript '0' stands for 9 = 0°. foj > 1 repre- 
sents field compression. Noj normalizes the field profile 
to 1. Po,/(a:) has an absolute value which is maximum 
in the center of the facet and zero at its edges. The con- 
stant ao./ has been introduced to account for idle region 
effects. 
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By decomposing H into components perpendicular and 
parallel to the ramp one realizes that Pqj < at least 
in the absence of idle regions. The screening currents in 
the idle regions reduce |Po,/(^)|- 

Similar as P90./, Bqj(x) can be regrouped as 



B 0J (x) = fi Q Hf Q j l + P 0J (x) 



(4) 



where Poj(x) = N j(Pqj(x) — (Poj)) has vanishing av- 
erage. 

If H is oriented within the (y,z) plane at an arbitrary 
angle 9 relative to the y axis the total flux density in the 
junction is 

B 9J {x) = B j(x) cos{9) + Bgo,f(x) sin(0) , (5) 

with B j(x) and B 90 j(x) as defined in Eqs. and @. 
The factors cos(0) and sin(0) arise from a decomposition 
of H into components parallel and perpendicular to the 
substrate plane. 

If the spatial dependences of Pqj(x) and Pgoj(x) are 
similar, the ratio p = \Pqj / Pqqj\ is about constant. 
Then, there is an angle 9h = — arctan(/o./p//go./) where 
the flux density Bgj penetrating the junction is homoge- 
neous along x and given by /jqH cos Ohfo,f(l—p/ho,f) ~ 
fj, H. 

The angle 9h might be referred to as "field applied par- 
allel to the ramp plane"—. In the absence of idle regions 
we expect Oh ~ r . In their presence 9h is reduced. We 
did not perform an explicit calculation, but a guess is 
to consider a field line, which starts at the YBCO ramp 

edge at y = [point (T) in Fig. 2] (a)] and touches the 
edge of the overlapping Nb film [point (2) in Fig. @] (a)], 
which, for <f> = 45° (the relevant angle for the multifacet 
system), is at the projected length y sa 4fim and z = 
125 nm. The corresponding angle is Oh = 1.8°. 

A similar argument will hold for a multifacet system. 
In the data shown below we determined the misalignment 
angle relative to 9h ("parallel" alignment, homogeneous 
field) as 6>* ff = 1.68° at <p = 45° and as 6>* ff = 0.85° at 
= 0°. We typically mount our samples with an offset 
angle 9 h < 1° relative to the substrate plane. Thus, we 
can put a limit Oh 2° which is fully compatible with 
the above estimate of 9h but rules out Oh = r . 

Finally, to account for the fact that Oh cannot be pre- 
cisely determined experimentally we introduce an angle 
0* = — Oh- By definition, the applied field is "parallel" 
for 9* = and "perpendicular" for 0* w = 90°. 

For the single facet case the flux density resulting from 
a field applied at an angle 9* is: 

B e ,j{x) = MoPcos(r) +(M 3 Hsm(0*)g X j(x), (6) 

with g±j(x) = /±,/[l + P±j(x)]. f±j denotes the field 

compression. As in Eq. ((2J, P±j-(x) « Pgo,/(x) is a spa- 
tially varying function with zero average, having its max- 
imum in the center of the junction plane at x = L/2. For 
the parallel field component (0* = 0° ) we have explicitly 



used that the flux density is homogeneous and essentially 
no flux compression occurs. 

We have not yet made use of the fact that field lines 
can penetrate the YBCO film freely in the junction area. 
Thus, in y direction the flux through the junction is not 
conserved but varies along y. Alternatively, assuming a 
constant flux density along y, this can be rephrased in 
terms of an effective junction thickness i e ff that varies 
along the y direction. Generally, t e g can be found via^ 



t e « — ^Au + Aybco tanh 



/ ^ybco \ 
V 2Aybco / 



+ A]Mb tanh 



( d 



Nb 



V2A 



Nb 



(7) 



Over the ramp area, the YBCO film thickness g?ybco 
grows along y from to ^ybco- Thus, t e g increases from 
~ 53 nm to ~ 85 nm over the ramp. For further calcula- 
tions we use a linearized Ansatz 



t e s(y) = icffO + At eff 



y_ 
w 



(8) 



with —0.5 < y/W < 0.5. The average effective thickness 
toffo is given by Eq. Q, using d Y BCO ~ ^ybco/ 2 - 

Ignoring self-field effects caused by the Joscphson cur- 
rent, the Josephson phase 7 is calculated via 



dj(x,y) 2tt 
ax 4>o 

The maximum supercurrent is obtained from 



(9) 



Ic(H) 



max 

7o 



dxdyj c {x)sin[j(x,y) -70] > (10) 



where integration has to be performed over the junc- 
tion area. The dependence of t c ff on y causes dephasing, 
which becomes severe when the flux difference (along y) 
over the junction width becomes on the order of $o/2- 
In the above estimate £ e ff differs by At e g/t e so ~ 0.5 from 
its average value i e fro (~69nm), and thus we expect the 
effect to become noticeable when the total flux through 
the junction becomes larger than about 4>o. However, 
the data discussed below for the 8 facet junction indicate 
At e s/t e ffo ~ 0.02. Also, using t c so as a free parameter, 
our data indicate a value for t e ffo which is close to 85 nm. 
I.e. only the "thick" part of the ramp seems to be rel- 
evant. Eq. (O assumed isotropic superconductors and, 
thus, the YBCO anisotropy could lead to modifications. 
Whether or not this solves the discrepancy is unclear to 
us. 

Let us summarize the central results of this section: 
We expect that in the junction plane the magnetic flux 
density along x follows the shape given by Eq. ((6]). B is 
homogeneous and equal to [IqH if the field H is applied 
under an angle Oh (~ 1° — 2° for our geometry) relative to 
the substrate plane, which we denote as "parallel field" 
(0* = — Oh = 0°). Perpendicular components of H lead 
to a compressed flux density which varies along x, having 
a maximum in the center of the junction at x = L/2. Due 
to the fact that flux lines can freely penetrate the YBCO 
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FIG. 5: (Color online), (a) Simulated in-plane magnetic field 
lines (black arrows) between two overlapping planar supercon- 
ducting layers (YBCO, red vertical lines, Nb, blue horizontal 
lines) separated by a distance <1l = 500 nm. The magnetic 
field was applied perpendicular to the layers. Inset shows field 
lines at expanded scales. For clearness blue lines are omitted, 
(b) Flux density profile Bgo(0 of the calculated in-plane mag- 
netic field B projected onto the zigzag line. The coordinate £ 
runs along the zigzag line at the edge of the YBCO film. 



film we expect that the total flux through the junction is 
not conserved but depends on y. This effect seems to be 
present in the zigzag junctions discussed below. However, 
its magnitude is at least an order of magnitude smaller 
than expected. 



2. Multifacet Ramp Junction 



plane and the field lines calculated in the plane parallel to 
the films and situated in the middle of the d^ = 500 nm 
gap between the films. The field was calculated along the 
edge of the YBCO film. B is strongest at the inner edges 
of the YBCO layer and weakest at the outer edges. A flux 
density profile Bgo(£,) — we use the index '90' to indicate 
that H has been applied perpendicular to the layers — of 
the calculated in-plane magnetic field B projected onto 
the zigzag line is shown in Fig.[5](b). £ is a curvilinear co- 
ordinate along the zigzag edge of YBCO. The projection 
is in units of HqH and the coordinate £ runs along the 
facets, — 0.5 < £/NL < 0.5. B 90 (t;) / [i H varies roughly 
sinusoidally, with one period per (0 - ir) - segment (by con- 
trast, one might have expected one period per facet, cf. 
Fig. @] (c)). £?9o(£) / HoH reaches a maximum value of 
about 4.5, which is in fact much less than the actual field 
compression (~ 100) found in experiment. It also turned 
out that the B component perpendicular to the layers is 
of the order of the in-plane component - a feature which 
we expect to disappear in a more realistic scenario. Fur- 
ther, simulating I C {H) with the profile of Fig. [5] (b) gave 
strong differences to the experimental interference pat- 
tern for flux values larger than $o per (0 - it) - segment. 

To simulate our zigzag junctions, as for the single facet 
case we use angles 9* = 9 — 9h and (f> to describe the 
angles between the applied field and the junction. In 
this case "perpendicular" refers to 9* — 90° and dif- 
fers from 9 = 90° by 9^. To account for this, below 
we use the index '_L' to denote perpendicular direction. 
For the flux density generated by the perpendicular com- 
ponent of the applied field we follow the two-facet peri- 
odic shape obtained by the 3D-MSLI simulations but, in 
order to be more flexible, generalize it in the following 
way. Each (0 - 7r) - segment is described by a field profile 
aj_j + Pj_j(x), with < x < 2L. We used 



( 7TX\ 

sin — 
V 2LJ 



(11) 



The resulting periodic pattern along £, i.e. the coordi- 
nate running along the zigzag line, is multiplied by an 
envelope function 



p±AZ) 



(12) 



extending smoothly across the whole junction. We 
parametrized Pj_ e via 



To get an idea of the field focused flux density pro- 
file inside the junction for H applied at 6 = 90°, we 
have simulated the flux density distribution around two 
overlapping planar superconducting layers located in the 
(x,y) plane being separated by a distance c?l along z, 
using 3D-MSLI^i, see hatched areas in Fig. [5] (a). The 
ramp was not included. The layer shape corresponded to 
the experimental situation and included N = 8 facets of 
length L — 10 /im. We used di = 500 nm. Smaller values 
of d^ led to convergence problems. The layer thicknesses 
were dyBCO = 65 nm and c?Nb = 100 nm, respectively. 
Fig. [S] (a) shows the geometry projected onto the (x, y) 



P 



±.e 



1 {nl) 



(13) 



The overall shape of B± (£) is given by 

£?x(0 = fx Hf ± N ± E(0 [a ±J + P ±if (£)] . (14) 

N± normalizes the spatial average of B±//j,oHf± to 1, 
f± denotes field compression. 

For arbitrary values of 9* and <f> the flux density B 
through facets oriented parallel to x is 



B e *(0 = B ± ($ s int 



Bh cos 9* cos 4> ■ (15) 
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The factors s'mO* and cos9* cosc/> arise from a projection 
of the applied field H. 

For facets oriented parallel to y one obtains 



«//(mT)-0.06 -0.03 0.00 0.03 



0.06 



B e , (£) = B ± (£) sin 9* + B« cos 9* sin < 



(16) 



with a homogeneous flux density B\\ — /j-qH along the 
zigzag line and B±(£) as defined in Eq. (fl"4|). Using this 
field profile (the actual parameters used are discussed in 
Sec. D, see also inset of Fig. [6] (a)) we solved Eqs. (|9|) and 
(fTtJ)) . with x replaced by £. The critical current density 
j c has been assumed to be constant in amplitude (ho- 
mogeneous junction). It changes sign between adjacent 
facets. 

Our model contains the 7 parameters otf, a e , aj_, e , 
At e g, t e ffo, /_L and 9h- Further, in experiment there is 
the offset angle 9 g relative to the substrate plane. As 
discussed in more detail in Sec. D, the shape of I C (H) for 
"perpendicular" fields is not very sensitive to small errors 
in 9 or 8* when 9h is small and thus, from this curve one 
obtains a/, a e , a± <e and At e g/t e go with good accuracy. 
From the field modulation period in perpendicular field 
one further obtains the product /_Lt G ffo- A priori, t e go is 
not known precisely. From the angle dependent I C (H) 
measurements one can identify the dead angle 0* d relative 
to 9h with high accuracy. However, 6* = (9 — Oh) = 0° is 
much harder to identify if, as in our case, I C (H) is mea- 
surable only over a small number of modulation periods. 
9^ depends on f±. Thus, f± needs to be determined by 
comparing a set of I C {H) curves near parallel orientation, 
finally also allowing to determine t a so and the offset angle 
9* s relative to 9h- 9h itself, however, is hard to determine 
and has an uncertainty which is on the order of the offset 
angle 9 a g. 



D. Results 
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FIG. 6: (Color online). Measured (red) and calculated (black, 
blue) interference patterns for (a) 8* = 90° and (b) <f> = 45° 
and 6* = 0° (calculation), 6* = —0.085° (measurement; offset 
angle 8* s = 1.684° subtracted). 8* — 9 — 9h = corresponds 
to "parallel" field orientation, as defined in the text. Model 
parameters are a e — 0.5, a±, e — 0.3, a/ = 0.7, a±j = 2.0, 
f± = 100, cf. Eqs. dTTJ), JTjj), {TJI and JTJJ. For the black 
curve At e f{/t e ffo = 0.02, for the dotted (blue) line At e g = 0. 
Ho/ f± — NLt c ffo/&o. Critical current in the calculated plots 
is normalized to I c o = b'cl-Aj, where j c is the critical current 
density and Aj is the junction area. Inset in (a) shows flux 
density profile B± normalized to (A)Hf±. 



Fig. HI (a) compares measured and calculated interfer- 
ence patterns for 9* = 90°. The horizontal scales of 
the calculated patterns are in units of Hq / f± , where 
Hq = f±NLt c ffQ/&o, with the field compression factor 
for perpendicular field components f± , the average effec- 
tive junction thickness t e so, the facet length L and the 
facet number N. The solid (black) line is for a relative 
variation of the effective junction thickness At e s/t e so = 
0.02, the dotted (blue) line is for At s = 0. The agree- 
ment between measured and calculated (for At e s /i e ffo = 
0.02) interference patterns is reasonable, although dif- 
ferences occur for normalized fields HJ±/Hq > 5. For 
Ai e ff/i e ffo > 0.02, at high fields the I c maxima are sup- 
pressed in comparison with the At e s/t e gQ = 0.02 case, 
while the I c minima are significantly lifted from zero, 
both in contrast to the measurements. Thus the mar- 
gin on Ai e flf/i e ffo is relatively narrow. For the (0-7r)- 
segments we used in Eq. (fTT|) aj = 0.7 and a±j- = 2.0. 
The power a/ was suggested by the 3D-MSLI simula- 
tions, cf. Sec. C. The calculated interference patterns, 
however did not depend strongly on this parameter. By 



contrast, the quite large value of a±j was necessary to 
achieve a reasonable agreement with experimental data. 
As a result, the periodic modulations of B±(£) are much 
shallower than suggested by the 3D-MSLI simulations. 
For the envelope function (Eqs. (fT2")) and (fl~3"|) ) we used 
a e = 0.5 and a±, e = 0.3. Normalization resulted in 7Vj_ 
= 0.34. The corresponding field profile is shown in the 
inset of Fig. [S] (a). Finally, by matching the abscissas 
of the theoretical and measured interference patterns we 
find t e sof± = 8.5 /an. A thickness f c gn = 69 nm calcu- 
lated from geometry (film thicknesses) corresponds to f± 
= 120. In fact, for the I C {H) data at different values of 
9* and <j) this value turned out to be somewhat too large. 
Best agreement was obtained for f± ~ 100, correspond- 
ing to t e ffo ~ 85 nm. 

We mention here that a periodic modulation of B±(£) 
with one period per facet rather than one period per 
(0 - tt) - segment also gave satisfactory agreement between 
calculated and measured interference patterns. Thus, 
from I C (H) we cannot unambiguously distinguish these 
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scenarios. 

Fig. [5] (b) shows measured and calculated interference 
patterns near parallel orientation for (j> = 45°. For the 
calculation (9* = 0°) we have assumed that the mag- 
netic flux density along the zigzag line of the multi- 
facet junction jBii (x) is homogeneous and not compressed, 
i.e. Bu(x) = fJ-oH. The field scale is given in units of 

ifo a/2- The factor v2 has been included to account for 
the fact that only a field [i^Hj 'y/2 is applied per facet. 
The main I c maximum in this normalization roughly 
appears at N/2. The dotted (blue) interference pat- 
tern is for Ai e ff = to emphasize the "ideal" case. 
Ates/t e so = 0.02 (black line) yields only minor devia- 
tions in the field range shown. The experimental curve 
is at an angle 9* = —0.085°. No curve has been mea- 
sured at 9* = 0.0° an angle, which is defined only after 
data analysis. Note that, in principle, the experimental 
and theoretical field scales arc linked for a fixed value of 
However, in order to compare the shape of I C (H), 
in Fig. [6] (b) we have additionally compressed the exper- 
imental field scale by a factor of 1.2 to account for the 
nonzero value of 9* . Further note that the experimental 
curve only covers a limited field range. This is the limit 
set by our experimental system. Comparing the curves 
at 9* = 90° and 6* = 0° one notices that for 9* = 0° the 
amplitudes of the secondary I c maxima situated between 
the two main I c maxima follow a U-shaped dependence 
while for 9* = 90° their amplitudes are about constant. 
Also, for 9* = 0° and absolute values of fields higher than 
the main I c maxima, the I c maxima are strongly reduced 
compared to 9* = 90°. These are the main differences in 
I C (H) at 9* = 0° and 9* = 90°. Thus, the overall dif- 
ference of the interference pattern at 9* = 90° compared 
to the ideal case is much less striking than in previous 
publications^—, presumably pointing to a much higher 
inhomogeneity of the critical current density in previous 
generations of ramp-zigzag - 7r J Js. 

Fig. [7] compares measured and calculated interference 
patterns for <p = 45° and -1.128° < 9* < 0.68°, where 
the angles quoted here and also below refer to the cal- 
culated patterns. In the experimental curves the offset 
angle 9* s = 1.684° relative to 9* = is not subtracted. 
The agreement is fair for all angles shown. The only addi- 
tional parameter required to fit the whole series of curves 
was 6** ff . It was determined by comparing the calculated 
interference pattern at the dead angle, cf. Fig. [7] (1), to 
the measured one, cf. Fig. [7] (k). For arbitrary </>, the 
dead angle (relative to 6* = 0) is given by 



<j) = 45° 



arctan 



2/i 



(17) 



with 9* w -0.4° for <j> = 45°. Near 9*, variations in 9* 
on the order of 0.01° already cause significant changes in 
I c vs. H so that, once f± is fixed, 6** ff can be determined 
very precisely. Note that the experimental I c vs. H 
pattern of Fig. [7J (k) slightly modulates around 2.5 /x A, 
while the calculated pattern in Fig. 0(1) is almost at zero 
current. This is presumably caused by residual fields in 
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FIG. 7: Experimental (left) and calculated (right) interference 
patterns for = 45° and -1.128° < 8* < 0.68°. 8* = 6-6 h = 
corresponds to "parallel" field orientation, as defined in the 
text. Offset angle 0* B = 1.684° relative to 8* — is not 
subtracted in the experimental graphs. Ho /f±~ NLt e m / <fro • 
Critical current in the calculated plots is normalized to I c q = 
[jcl-Aj, where j c is the critical current density and Aj is the 
junction area. Model parameters are a e = 0.5, a± <e — 0.3, 
a f = 0.7, a X j = 2.0, f± = 100 and At cH /t cHa = 0.02, cf. 
Eqs. <[TTJ), (T2J>, CG2> and ([14]). 



the cryostat, causing also the shift in the I c minimum 
relative to H = in Fig. [7](m). Also, the 9* dependence 
of the interference patterns strongly depends on f± . This 
parameter cannot be altered by more than some 5 % from 
100 without substantial degradation of the fit qualitiy. 

To further confirm that the profile -B(£) - plus the 
assumption of homogeneous junction parameters - de- 
scribes the situation well we have also varied (j>. Fig. [8] 
compares data and calculations for = 0° and var- 
ious values of 9* around 9* = 0°. For the case of 
4> = 9* = 0° calculations yield critical current main max- 
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0= 0° 
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FIG. 8: Experimental (left) and calculated (right) interference 
patterns for = 0° and -0.846° < 9* < 0.823°. 9* = 6-6 h = 
corresponds to "parallel" field orientation, as defined in the 
text. Offset angle 8* s = 0.846° relative to 9* — is not 
subtracted in the experimental graphs. Ho/f± = NLt c m)/&o- 
Critical current in the calculated plots is normalized to I c q — 
\j c \Aj, where j c is the critical current density and Aj is the 
junction area. Model parameters are a e = 0.5, a± t£ — 0.3, 
a ; = 0.7, a±j = 2.0, f± = 100 and At eB /t eB0 = 0.02, cf. 
Eqs. QT), dl, (O and pi). 



ima I c = 0.5/ c o whenever the flux through the facets 
oriented along x equals a multiple of <&o- Then the crit- 
ical current of these facets cancel and I c is given by the 
sum of the critical currents of the facets oriented along y. 
Curves of Fig. [5] (a)— (f) show the appearance of the first 
and second main I c maximum, which grows in amplitude 
for 9* — » 0°. Although the maximum field provided by 
our setup was not sufficient to resolve the main max- 
ima for 6* = 0° the evolution of the first and second 
I c main peaks are clearly visible for 6* < 0.336°, with 
a good agreement between experimental and theoretical 



patterns. By lowering 6* one again runs through a dead 
angle (6*2 = —0.253°) and obtains somewhat strangely 
looking interference patterns for —1° < 6* < 6* d . 

Finally, we address the effect of self fields generated by 
the Josephson currents. Such effects become prominent 
when the junction is 2-3 times longer than the Joseph- 
son length A j. To estimate Aj we first ignore idle re- 
gions and assume that the supercurrent flows homoge- 
neously across the junction area Aj — NLW « 37/im 2 . 
We then find j c0 = I c0 /Aj w 140 A/cm 2 . Using Xjo = 
[$o/27TMoJcrfoff] ' 5 , with the effective magnetic junction 
thickness^ 

d e s = t Au + Aybco coth ( ^ YBCO J -f A Nb coth [ ) 

V-^YBCO/ \ANb/ 

w 470 nm (18) 

we find A,/o ~ 20 (im and the normalized junction length 
I = NL/Xj sa 4. The idle region effect^ 2 .^ leads to 
an increased Josephson length, Aj ; = S ■ Ajo, with 
S = (1 + (d eS /d eS! i)(Wi/W)) - 5 . Wi and 4ff,i, respec- 
tively, are the width and effective magnetic thickness of 
the idle region. With d c g t \ = 535 nm and Wi = 3 fjja 
one obtains 5 f=a 2.6 and Xj^ w 50 /iin. Thus, / 1.6, 
justifying the short junction approach taken above. An 
ambiguity, however, arises from the problem to refer the 
measured I c to the proper junction area. Above, we have 
used the whole ramp area. Alternatively, assuming that 
the current is dominantly carried by the in-plane cur- 
rents on the YBCO side one might refer to a junction 
area which is projected perpendicular to the substrate 
plane. Then, j c increases by a factor of 7 and Aj de- 
creases by a factor of 3, bringing the junction closer to 
the long junction regime. To distinguish these scenarios 
we simulated I C (H) based on the sine-Gordon equations 
for the case of 6* = 90° and various values of the nor- 
malized junction length /. For the flux density profile the 
same shape as in the inset of Fig. [5] (a) has been used. 
Fig. [9] shows simulations of I C (H) patterns for different / 
and the calculation using the short junction model. Devi- 
ations from the short junction model that are incompat- 
ible with our experimental data become prominent near 
the main I c maxima for I > 3 (the case of I = 4 is shown 
in the graph). By contrast, I C {H) for / = 1.6 is almost 
indistinguishable from the short junction model and in 
agreement with the above estimate, using the ramp area 
instead of its projection. 



E. Conclusions 

In summary, we have investigated the interference pat- 
terns I c vs. H for a 8 facet YBCO-Au-Nb zigzag-ramp 
Josephson junction, with a facet length of 10 /jm. The 
angle 9 between the substrate plane and the applied field 
H was systematically varied and we also discussed two 
in-plane angles <j) between H and the facet orientation 
(cf) = 0° and (f> = 45°). All interference patterns could 
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FIG. 9: (Color online). I c vs. H (6* = 90°) patterns calcu- 
lated with the short junction model for At e g = (solid line) 
and simulated for different lengths / (circles) using StkJJ^i 
in comparison with measurement (dotted line). Other model 
parameters are the same as in Figs. [6] [7] and [8] Deviations 
between the calculated and the simulated (for I = 4) pattern 
are visible in the zoomed inset. 



be understood from the fact that a field component per- 
pendicular to the substrate causes a strong and spatially 
varying contribution to the flux density profile inside the 
junction. The junction by itself - admittedly our best 
junction - seems to be very homogeneous, with an essen- 
tially constant critical current density inside each facet. 
Particularly, no j c asymmetry between facets oriented 
along the in-plane x and y directions were observable. 
We have expected a strong dephasing effect on I c vs. H 
due to the fact that the flux penetrating the junction is 
not preserved. The effect is present but at least an order 
of magnitude weaker than expected. Our investigations 
also showed that H is applied "parallel" to the junction 
(in terms of an homogeneous flux density profile) for an 
angle 9h ^ 2° and not for a field orientation about par- 
allel to the ramp angle 9 r = 8°. These findings may 
contribute new knowledge to the general physics of ramp 
junctions. 



Further, due to strong compression of the perpendic- 
ular field component by about a factor of 100, only for 
angles \9*\ = \0 — 9h\ 0.1° the interference patterns 
were "ideal" in the sense that the flux density in the 
junction is essentially homogeneous. The effective junc- 
tion thickness is only about 85 nm, leading to enormous 
fields (~ 35 mT) that are required to produce a flux quan- 
tum per (0 - 7r) - segment for the case of 9* = 0° . We also 
demonstrated, that there is a dead angle 9* d very close to 
"parallel" field orientation, where tremendous changes in 
I C (H) occur. These properties make it extremely hard to 
study zigzag junctions in "parallel" field configuration. 
Due to idle regions the Josephson length of our junction 
was on the order of 50 /jm, requiring junction lengths of 
several 100 /im to study long junction effects. For such 
junctions the field compression will be even bigger than 
in our case, and, thus, the field alignment must be much 
better than 0.1° relative to 9* = to achieve a homo- 
geneous flux density. Thus it seems that, for a study of 



(semi)fluxon physics 



3.6.7.15 



or (^-junction effect; 



,18-21 



zigzag junctions magnetic fields should be oriented per- 
pendicular to the substrate plane and the corresponding 
flux density profiles should be taken into account rather 
than being avoided by parallel field alignment. 

We finally note that dead angles, as described here also 
appear in other types of Josephson junctions or Joseph- 
son junction array o 30 ' 32 ' 35 . Whenever there is a strong 
compression of the perpendicular field component this 
angle may be very close to the parallel field orientation, 
making measurements in "parallel" field unreliable, fa- 
voring the perpendicular field orientation. 
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